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GAUSSIAN DISTRIBUTION OF THE NUMBER OF SUMMANDS IN GENERALIZED 
ZECKENDORF DECOMPOSITIONS IN SMALL INTERVALS 


ANDREW BEST, PATRICK DYNES, XIXI EDELSBRUNNER, BRIAN MCDONALD, STEVEN J. MILLER, 

KIMSY TOR, CAROLINE TURNAGE-BUTTERBAUGH, AND MADELEINE WEINSTEIN 

ABSTRACT. Zeckendorf’s theorem states that every positive integer can be written uniquely as a sum 
of non-consecutive Fibonacci numbers F n , with initial terms F\ = 1.7*2 = 2. Previous work proved 
that as n —> 00 the distribution of the number of summands in the Zeckendorf decompositions of 
to £ [F n ,F n+ 1 ), appropriately normalized, converges to the standard normal. The proofs crucially 
used the fact that all integers in [F n ,F n+ 1 ) share the same potential summands and hold for more 
general positive linear recurrence sequences {G n }. 

We generalize these results to subintervals of \G n . G n + 1 ) as n —> 00 for certain sequences. The 
analysis is significantly more involved here as different integers have different sets of potential sum¬ 
mands. Explicitly, fix an integer sequence a{n) —> 00 . As n —>■ 00 , for almost all to £ [G n ,G n +i) 
the distribution of the number of summands in the generalized Zeckendorf decompositions of integers 
in the subintervals [to,to + G ct ( I j)), appropriately normalized, converges to the standard normal. The 
proof follows by showing that, with probability tending to 1, to has at least one appropriately located 
large gap between indices in its decomposition. We then use a correspondence between this interval 
and [0, G, l:(n) ) to obtain the result, since the summands are known to have Gaussian behavior in the 
latter interval. 
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1. Introduction 


1.1. Background. Let {F n } denote the Fibonacci numbers, normalized so that = 1, F 2 — 2, and 
F n+ 1 = F n + F n _ i. One of the more interesting, equivalent definitions of the Fibonacci numbers is 
that they are the unique sequence of positive integers such that every positive number has a unique 
legal decomposition as a sum of non-adjacent terms!] This equivalence is known as Zeckendorf’s 
theorem (Zell . Once we know a decomposition exists, a natural question to ask is how the number of 
summands varies. The first result along these lines is due to Lekkerkerker HLekH . who proved that the 
average number of summands needed in the Zeckendorf decomposition of an integer m e | F n . F n+l ) 
is -q-j- + 0(1), where p = 1+ 2 ^ . the golden mean, is the largest root of the Fibonacci recurrence. 
More is true, and many authors have shown that the distribution of summands of m € | F n . F n+] ) 
converges to a Gaussian. These results have been extended to a variety of other sequences. There 
are several different methods of proof, from continued fractions to combinatorial perspectives to 
Markov processes; see HB-AM1 |Day| IDDKMV1 iDGl IFGNPT1 IGH IGTNP1 iKel IKKMW1 ICTl iLenl 
IMWll IMW21 IStell ISte2H for a sampling of results and methods along these lines, (All ICHMNll 
ICHMN21ICHMN31 IDDKMMVl IDDKMVH for generalizations to other types of representations, and 
(BBGILMTl IBTLMPl for related questions on the distribution of gaps between summands. 

The analysis in much of the previous work was carried out for m e j F n . F n+] ) (or, for more 
general sequences {G n }, for m e [G ni G n+ 1 )). The advantage of such a localization] is that each 
m has the same candidate set of summands and is of roughly the same size. The purpose of this 
work is to explore some of the above questions on a significantly smaller scale and determine when 
and how often we obtain Gaussian behavior. Note that we cannot expect such behavior to hold for 
all sub-intervals of [F n , F n+ 1 ), even if we require the size to grow with n. To see this, consider the 
interval 

\F 2 r 1 + F n + F n _ 2 + • • • + Fy n i/4j , +2n + F n+ 1 + Fy n i/4j). (LI) 

The integers in the above interval that are less than F2n + F n+ 1 have on the order of n/2 summands, 
while those that are larger have at most on the order of n 1/4 summands. Thus the behavior cannot be 
GaussianH 

In f BDEMMTTWl . we proved Gaussian behavior for the number of summands in the Zeckendorf 
decomposition for almost all small subintervals of [F ni F n+1 ); in this work we generalize to other 
sequences. Henceforth {G n } will denote a positive linear recurrence sequence: 

Definition 1.1. A sequence {G n }f =l of positive integers is a Positive Linear Recurrence Sequence 
(PLRS) if the following properties hold: 

(1) Recurrence relation: There are non-negative integers cl such that 


G n +i 

with L, ci and c L positive. 


— c\G n + • • • + CiG„ + i_i, 


(2) Initial conditions: G\ = 1, and for 1 < n < L we have 

G n +1 = ciG n + C 2 G n -i + • • • + c n G\ + 1. 


( 1 . 2 ) 


(1.3) 


'The requirement of uniqueness of decomposition forces us to start the sequence this way. 

“As the sequence { F n } is exponentially growing, it is easy to pass from m in this interval to m G [0, F n ). 
Though in this situation it would be interesting to investigate separately the behavior on both sides. 
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A decomposition Y^T=\ a iG m +i-i of a positive integer N (and the sequence {«,:}(= if is legal if 
a\ > 0, the other a* > 0, and one of the following two conditions holds: 

• Condition 1: We have m < L and ai = Cifor 1 < i < m. 

• Condition 2: There exists s G {0, .... L} such that 

ai = ci, a 2 = c 2 , , GLj-i = c s _! and a s < c s , (1.4) 

a s+ 1 , ..., a s+ c = 0 for some i > 0, mzd {biYfLi (with hi = a s+ e +i ) is legal. 

If a iGm.+ \-i is a legal decomposition of N, we define the number of summands (of this 
decomposition ofN) to be a 4 + • • • + a m . 

Informally, a legal decomposition is one where we cannot use the recurrence relation to replace 
a linear combination of summands with another summand, and the coefficient of each summand is 
appropriately bounded. For example, if G n+ 1 = 2 G n + 3G n _i + G n _ 2 , then G 5 + 2 G 4 + 3G 3 + G 4 is 
legal, while G 5 + 2G 4 + 3G 3 + G 2 is not (we can replace 2G 4 + 3G 3 + G 2 with G 5 ), nor is 7G 5 + 2 G 2 
(as the coefficient of G 5 is too large). Note the Fibonacci numbers correspond to the special case of 
L = 2 and c 4 = c 2 = 1. 

Earlier work (see, for example, I MWll IMW21 ) proved that if {&'„} is a PLRS then we again 
have unique legal decompositions and Gaussian behavior for the number of summands from m G 

[G n , G n+ 1 ). 

1.2. Main Result. 

Fix an increasing positive integer valued function a(n) with 

lim a(n ) = lim (n — a(n)) = 00 . (1.5) 

n—>00 n —^00 

Our main result, given in the following theorem, extends the Gaussian behavior of the number 
of summands in Zeckendorf decompositions to smaller intervals. Note that requiring m to be in 
[G n , G n+ i) is not a significant restriction because given any m, there is always an n such that this 
holds. 

Theorem 1.2 (Gaussianity on small intervals). Let {G n } be a positive linear recurrence sequence 
with recurrence 


G n + 1 — ciG„ + c 2 G n _i + ■ ■ ■ + clG u+ \-l, (1.6) 

where we additionally assume c\ > c 2 > ■ ■ ■ > cl > 1. Choose a(n) satisfying (fT3T) . The 
distribution of the number of summands in the decompositions of integers in the interval \m. rn + 
G a ( n )) converges to a Gaussian distribution when appropriately normalized for almost all m G 

[G n , G n+ 1 ). 

In C2 we derive some useful properties of Zeckendorf decompositions, which we use in CH to 
prove Theorem ll.21 The reason for our extra condition on the cf s surfaces in Lemma [2T4l where this 
constraint forces a truncated legal decomposition to remain legal. We conjecture that Theorem 11.21 
holds more generally. 
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2. Preliminaries 


In order to prove Theorem 11.21 we establish a correspondence between the decompositions of 
integers in the interval [m,m + G a { n )) and those in [ 0 . G n ( n) ). Throughout this paper, when we 
write an interval [a, b], we mean the integers in this interval. We first introduce some notation. Fix a 
non-decreasing positive function q(n), taking on even values, such that 

q(n) < n — a(n), q(n ) = o ( \J a(n )) , lim q(n) = oo; ( 2 . 1 ) 

the reason for the second condition is to allow us to appeal to known convergence results for a related 
system, as q{n) will be significantly less than the standard deviation in that setting. 

For m G [G n , G n+ 1 ) with decomposition 


define 



( 2 . 2 ) 


C'r(m) 

C 2 (m) 

C 3 (m) 

Note that each a t G {0,1} for all 1 
decomposition of m. That is, let 


• (Oi, O 2 j •••) ®a(n))i 

• (®a(n)+l) ®a(n)+g(n)) > 

: = (Oa(ra)+g(n)+l) ■■■, On)- (2.3) 

< i < n. Let s(m) be the number of summands in the 


s(m) := (2.4) 

3 = 1 


Similarly, let Si(m), S 2 (m), and s 3 (m) be the number of summands contributed by Ci(m), C 2 (m), 
and C 3 (m ) respectively. 


Lemma 2.1. Let x G [m, m + G a ( n )). If there are at least 3 L consecutive 0’s in (^(m), then C 3 {x) 
is constant, and hence s 3 (x) is constant as well. 


Proof. Assume there are at least 3L consecutive 0’s in C 2 {rn). Then for some k e [a(n) + 3L, a(n) + 
q(n)), we have afc_ 3 L+i = a k - 3 L+2 = • • • = a k = 0. Let m' denote the integer obtained by truncating 
the decomposition of m at a k - 3 LGk- 3 L■ Then ml < G k - 3 L+i- Since G a ( n ) < G k ~ 3 L, it follows that 
for any h < G a ( n ) we have 

m' + h < G k - 3 L+i + G k - 3 L 

< Ci G k -3L+1 + c 2 G k - 3 l + ■ • • + c^Gfc -2L+2 ~ G k - 3 L+2 5 (2.5) 

and thus the decomposition of vn! + h has largest summand no greater than G fc _ 3L+2 . Then since 
3L — 2 > L, the Zeckendorf decomposition of rn + h is obtained simply by concatenating the 
decompositions for m — m' and vn! + h. Hence C 3 (m + h) = C 3 (m — m') = C 3 {m). □ 


With this lemma, we see that the distribution of the number of summands involved in the decom¬ 
position of x G [m, m + G a (n)) depends (up to a shift) only on what happens in C\(x) and C 2 (x), 
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provided that there is a gap between summands of length at least 3 L somewhere in C 2 (m ). In light 
of this stipulation, we will show the following items in order to prove our main theorem. 

• With high probability, m is of the desired form (i.e., there is a gap between summands of 
length at least 3 L in C 2 (m)). 

• When m is of the desired form, the distribution of the number of summands involved in 
C±(x ) for x G [m, m + G a ( n )) converges to Gaussian when appropriately normalized. 

• The summands involved in C 2 {x) produce a negligible error term (i.e., there are significantly 
fewer summands from C 2 (x) than there are from C\ (rn)). 

We address the first point with the following lemma. 

Lemma 2.2. As n —> oo, with probability 1 + o(l) there are at least 3 L consecutive 0’s in C 2 (m) if 
m is chosen uniformly at random from the integers in \G n . G n+ 1 ). 

In HBILMT1 . Bower, Insoft, Li, Miller and Tosteson analyze the largest gap between summands in 
legal decompositions, assuming only that the c^’s are non-negative (a less restrictive condition than 
we have). They prove that for m G \G n , G n+ 1 ) as n —> oo with probability 1 + o(l) there is a gap 
of length C logn for some constant C > 0. This almost proves what we need, the only difficulty is 
we need the gap to be in C 2 (m). Instead of modifying their technical approach, in Appendix [A] we 
derive a difference equation which proves as n —>• oo with probability 1 we have a gap of any fixed 
finite size; while we prove that such a gap exists somewhere, the argument there is trivially modified 
to apply to just C 2 {m) (as if we have no summands in that region then the result is trivially true!), 
and completes the proof of the above lemma. 

Until this point we have not used the condition that c\ > c 2 > ■ ■ ■ > cl- The following lemma 
illustrates the necessity of this condition. 

Lemma 2.3. If c\> c 2 > ■ ■ ■ > cl, then 

n n 

^ cijG n+1 ^j is legal a ? G n+w is legal. ( 2 . 6 ) 

i =i j =2 

Proof If the first sum is legal, then for some s we have a\ — Ci, a 2 = c 2 , ■ ■ ■ : a s < c s . Then we 
necessarily have a 2 < ci, a 3 < c 2 , • • • a! s < c s /„i for some s' < s. We can repeat this process, 
noting that a block ends at a s . Therefore, the decomposition is legal. □ 

Corollary 2.4. If c\ > c 2 > • • • > c r , then legal decompositions are chosen by the greedy algorithm. 

Proof Suppose not, so that YTj=i a jG n +i-j is a legal decomposition with )T'' ; , + | a 3 G n+ \_j > 
G n+ i_fc. By an iterated application of the previous lemma, the above sum is a legal decomposition, 
and it is known that legal decompositions of integers in [G n , G n+ 1 ) all have largest term G n , which 
gives a contradiction. □ 

Assuming m is of the desired form (meaning there is a gap of length at least 3 L in its decomposi¬ 
tion), we now consider the distribution of s(x) for 16 [m,m + G Q ( n )). 

Lemma 2.5. If m has at least 3L consecutive 0’s in C 2 (m), thenforallx G [m,m +G a ( n )), we have 

0 < s(x) — s 3 (m) — s(t(x)) < Kq(n ), (2.7) 
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where K := 


given by 


maxj Cj, and t(x) denotes the bijection 

t : [m, m + G a ( n )) —> [0, G a ( n )) 


t(m + h ) 


m 0 + h, 

m 0 + h — G a (n ), 


if m 0 + h < G a (n ), 
if m 0 + h > G a ( n ), 


where m 0 is the sum of the terms in the decomposition ofm truncated at a a ( n )_iG Q („)_i. 


( 2 . 8 ) 


(2.9) 


Proof. First, note that the number of summands in the decomposition of x with indices i G [o(n), 
a(n) + q(n)) must be less than Kq(n). Next, we verify that t from above in fact is a bijection, for 
which it suffices to show that it is an injection. Note that if t(m + hf) = t(m + h 2 ), then since 
\hi — h 2 \ < G Q (n)> we must have m 0 + hi = m 0 + h 2 , so we may conclude that t is injective, hence 
bijective since the domain and codomain are finite. For any x G [m, m + G n ( n] ) , the decompositions 
of t(x) and x agree for the terms with index less than a(n) by virtue of Corollary 12.41 Furthermore, 
the decompositions of x and m agree for terms with index greater than a(n) + q(n). Therefore, the 
number of summands in the decomposition of x with indices i G [a(n), a(n) + q(n)) is equal to 
s(x) — s 3 (m ) — s(t(x)). Combining this with our initial observation, the lemma now follows. □ 


As a result of this lemma, the distribution of s(x) over the integers in [m, m + G a ) is a shift of 
its distribution over [0, G a ^), up to an error bounded by q(n). With this fact, we are now ready to 
prove the main theorem. 


3. Proof of Theorem II .21 

We now prove our main result. We assume below that {G n } is a PLRS with ci > c 2 > • • • > cl > 
1, that a(n) obeys (11.51) . and that q(n) obeys (12.11) . 

Proof of Theorem 17721 For a fixed m G [G ni G n+ 1 ) with at least 3 L consecutive 0’s somewhere in 
C 2 (m ) (which is true for almost all m as n -x oo), we define random variables X n and Y n by 

:= siHf), Y n := s(H 2 ), (3.1) 

where Hi is chosen uniformly at random from [m, m + ) and II 2 = Hi— m (and thus is chosen 

uniformly at random from [0, G a ( n ))). Let 

K ■= -^ n -(X n -E[X n ]), (3.2) 

o x {n) 

and 

K ■■= - E[XJ), (3.3) 

where a x (n) and a y (n) are the standard deviations of X n and Y n , respectively, so that X n and Y n 
are normalized with mean 0 and variance 1. It is known that Y' n converges to a Gaussian distribution 
with mean and variance of order a(n) (see, for example, IMW1] ), and we claim that X' n converges 
to a Gaussian distribution as well. 
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Note that X n = Y n + C + e(n), where C is a constant and e(n) is an error term with \e(n)\ < 
Kq(n). From the assumption that q(n) = o we have q(n)/a y (n ) —> 0. Therefore, we 

have 


Var[X n ] = Var[F n ] + Var[e(n)] + 2Cov[y n , e(n)] 

< Var[y n ] + AK 2 q(n) 2 + 2E[(Y n - E[y n ])(e(n) - E[e(n)])] 

< Var[y n ] + AK 2 q(n) 2 + 2E[|y n - E[y n ]| • |e(n) - E[e(n)]|] 

< Var[y n ] + 4 K 2 q(n) 2 + AKq(n)E[\Y n - E[F n ]|] 

< Var[y n ] + AK 2 q(n) 2 + AKq(n)E[{Y n - E[F n ]) 2 ] 1/2 
= Var[y n ] + AK 2 q(n) 2 + AKq{n) + AKq(n)a y {n) 

= Var[y n ](l + o(l)), (3.4) 


Hence a x (n) = a y (n)(l + o(l)). Note that \X n — Y n — C\ < Kq(n), so \X' n a x (n) — Y^(T y {n)\ = 
0(q(n)), and thus 


ix-y'i = 


x , °~x(w) _ , 

^n ( \ 1 n 

a v (n) 


(l + o(l)) = O 


a y (n) 


= o(l). 


(3.5) 


Let A n and B n be the cumulative distribution functions for X' n and Y' n , respectively. By Lemma 
12.51 combined with the above bound on \X' n — Y^\, we have 


B n (x-o( 1)) < A n (x) < B n (l + o(l)). 


(3.6) 


Since {B n } n converges pointwise to the cumulative distribution function for a Gaussian distribu¬ 
tion, say B(x), and it is easy to show that if /3 n —> 0, then B n (x + (3 n ) —>■ B(x) since B is continuous 
and B n is monotone non-decreasing, it follows that {A n } n also converges pointwise to B(x). This 
completes the proof. □ 


4. Conclusion and Future Work 

By finding a correspondence between generalized Zeckendorf decompositions in the interval 
[m, m + G a (n)) and in the interval [0, G a ( n )), we are able to prove convergence to Gaussian behavior 
on many sub-intervals. The key step is to show that almost surely an integer m chosen uniformly 
at random from [G n , G n+ 1 ) permits the construction of a bijection onto the interval [0, G a ( n )). Our 
results then follow from previous work on the Gaussian behavior of the number of generalized Zeck¬ 
endorf summands in this interval. 

In the future, we plan to extend our results to more general sequences. The first natural candidate 
is to remove the assumption on the c t ’s among the positive linear recurrence sequences we study. 
Other interesting topics include the signed decompositions (or far difference representations) where 
both positive and negative summands are allowed (see ITATI IDDKMVI ), the /-decompositions of 
IlDDKMMVl . and some other recurrences where the leading coefficient is zero (which in some cases 
leads to a loss of unique decompositions), such as UCHMN1 1 1CHMN21 Cl 1MN3I . 
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Appendix A. Elementary Proof of Moderate Gaps 

A crucial ingredient in our proof is that almost surely the Zeckendorf decomposition of an m E 
[G n , G n+ 1 ) has a gap of length Z or more for some fixed Z > 3 L. This follows immediately from the 
work of Beckwith, Bower, Gaudet, Insoft, Li, Miller and Tosteson [BBGILMTJ (see also HB-AMlI l. 
who showed that almost surely the longest gap is of the order log n. It is possible to elementarily 
prove this result by deriving a recurrence relation for these probabilities and analyzing its growth 
rate directly, which we do below both in the hopes that it might be of use in related problems, and 
also to keep the argument elementary. 

Theorem A.l. Let { G n } be a PLRS with recurrence relation with positive integer coefficients: 

G n +i = c\G n + • • • + clG u _l+i, Ci, L > 1 . (A.l) 

Fix an integer Z > L and let H n+1 be the number of integers m E [ 0, G n+ f) such that m’s legal de¬ 
composition does not have a gap between summands of length Z or greateru Then H n /G n converges 
to zero exponentially fast. 

Proof We first find a recurrence relation for {H n }. Consider n much larger than L + Z, and let 
m E [0, G n+ 1 ) be arbitrary. We count how many m have a legal decomposition with no gap between 
summands of length Z or more by looking at the possible beginning strings of m’s decomposition. 
Specifically, we look at how often G n , G n _i,..., G u -l+i occur. As the analysis is trivial when 
L = 1 we assume L > 2 below. 

• Case n: At most c\ — 1 summands of G n : We may have 0, 1,..., c\ — 1 or c\ occurrences 
of G n . If we have zero Gf s, then the number of possible completions of the decomposition 
such that all gaps are less than Z is, by definition, H n . If we have 1, ..., C\ — 1 summands 
of G n then the number of completions is H n — H n _ z (we must subtract H n _ z as we have a 
summand, and we cannot have Z or more non-chosen summands now). Thus the contribution 
from this case to H n+ 1 is 

H n + (ci — l)(H n — H n _ z ) = Ci(H n — H n _ z ) + H n _ z . (A.2) 

• Case ra—l: At most c 2 — 1 summands of G n _i: To be in this case, we first have C\ summands 
of G n . If we have no G n _i terms in our decomposition, then we have H n _i — H n _ 1 _^ z _ i) = 
H n _ i — H n _ z ways to complete the decomposition (as we have zero summands, we must be 
careful and avoid not taking any of the next Z — 1 summands). If we have 1,..., c 2 — 1 copies 
of G n -1 then there are ] — II„-\- Z ways to complete the decomposition as desired. Thus 
this case contributes to H n+ 1 

(-ffn-l — H n _ z ) + (c 2 — 1) [H n —\ — H n _i_ Z ) 

= c 2 (H n _ 1 — H n _i_ z ) + ( H n _i_ z — H n _ z ). (A.3) 

• Case n—£: At most q— 1 summands of G„_/; + 1 : Similar to the earlier cases, if we have zero 

copies of G n -e + i then there are H n _£ — = H n _i — H n _i +1 _ z ways to complete 

4 Note that we only care about a gap of length at least Z between adjacent summands; thus our decomposition may 
miss many elements before the first chosen summand. It is also easier to first count in [0, G n _|_i) and later shift to 
[G ni G n + 1 ) by subtraction. 
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the decomposition, while if there are 1,..., eg — 1 copies of G n _g + \ there are H n _g — H n _g_ z 
possibilities. Thus the total contribution to H n+1 is 

(■ H n ~e — H n _g + i_z) + (eg — 1 ){H n -g — H n _g_z) 

= ce(H n -f — H n _£_z) + (H n _e_z ~ -ffn-^+i-z)- (A.4) 


If i < L we continue to the next case with now eg copies of G n -g + 1 , while if £ = L then the 
process terminates, as the recurrence relation and our definition of legality means we would 

replace this beginning string with CiG n H-b c L G n _ L+ 1 with G n+ 1 , contradicting both the 

fact that we have a legal decomposition and that our number is less than G n+1 . 

Combining the above, we obtain the recurrence for {H n } (notice that we have a telescoping sum) 


C\(H n H n _ z ) + • • • + cz{^H n _z-s r \ H n _i / J r \_z^) -t- H n —L+i—z• (A.5) 


Setting H n+ i = H n+1 — H n , we see H n+ 1 counts the number of integers in [G n , G n+ \) with no gaps 
of length Z or more, and using (IA.5I) with n and n — 1 gives us (IA.5I) with each Hk replaced with H^. 
As {H n } is a strictly increasing sequence, the largest eigenvalue of its recurrence must be greater 
than 1. From the Generalized Binet Formula expansion for H n (see, for example, Theorem A.l of 
[BBGILMT]), there is some constant a G (0,1) such that H n _g_z > <y.H n _g. Letting 


we find 


eg 


cg( 1 — a) if 1 <i<L 
c L if £ = L, 


H n + i < C\Hx + ■ ■ ■ + clH u _i, +1 . 


(A. 6 ) 


(A.7) 


Notice this is almost the same recurrence as that of G n , the only difference being that at least one 
of the coefficients q is smaller than eg. To prove { II,,} grows exponentially slower than {G n }, we 
instead study the sequence {H n } which satisfies the recurrence 

H n +1 = c\H n + • • • + clH u _l+ i, (A.8) 

as H n < H n . 

We use many standard properties of the Generalized Binet Formula expansions below; see Theo¬ 
rem A. 1 of | BBGH.MTil for statements and proofs. The solution to a linear recurrence of fixed, finite 
length is of the form 

+ • • • + Pi, o-^i + • • • + fik,r k n rk ^k + ''' + Pk,o^k- (A.9) 

As the Ci s are positive, there is a unique positive root A > 1 for {G n }, and all other roots are less 
than 1 in absolute value. Thus G n = (3\ n + 0(|n r A 2 | n ) for some 0 < |A 2 | < A (if (3 = 0 then G n is 
exponentially decaying). As all the coefficients of the recurrence for {H n } are positive, there is also 
a unique root of largest absolute value. 

We claim to < A. Clearly to < A as this sequence grows slower; if they were equal then to would 
be a root of both characteristic polynomials, and we would find 

oj — C\0J T ■ ■ ■ T Cl < C\iu T ■ ■ ■ T cl = ^ ■ 

a contradiction. Thus H n grows exponentially slower than G n , which completes the proof. 


(A. 10) 
□ 
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